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We analyze the direct detection rate of minimal supersymmetric neutralino dark matter in germanium, sapphire 
and sodium iodide detectors, imposing cosmological and accelerator bounds including those from b — > S7 decay. 
In contrast with several other recent analyses we find models with light charged higgsinos and large stop mixing 
in which the counting rate in solid state detectors exceeds 10 events/kg/day. 
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The recent observation of the b — > S7 decay 
has stirred interest in the possible bounds obtain- 
able for supersymmetric models that contribute 
to this process Some authors |^,|| ana- 

lyzed the consequences of this and other accel- 
erator bounds on the predicted rates in experi- 
ments searching for neutralino dark matter in the 
galactic halo. They claimed that allowed rates 
are small, largely irrelevant for present-day dark 
matter detectors. 

We performed a detailed study of the allowed 
supersymmetric parameter space. In contrast 
with , we found models in which the integrated 
counting rates are not at all small but as high as 
to have already been probed (and excluded) by 
current negative dark matter searches. Here we 
present an overview of our approach and results, 
and refer the reader to for further details. 

We work in the framework of the minimal 

= 1 supersymmetric extension of the standard 
model (MSSM) We include one-loop radiative 
corrections in the Higgs sector according to the 
effective potential approach. For simplicity, we 
make a simple ansatz of a universal (weak-scale) 
sfermion mass parameter m j. This ansatz implies 
the absence of tree-level flavor changing neutral 
currents in all sectors of the model. Our remain- 
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ing arbitrary parameters are the soft supersym- 
metry breaking trilinear couplings Af, and At , the 
ratio of the two Higgs vacuum expectation values 
tan/3, the pseudoscalar mass myi, the higgs (ino) 
mass parameter /i and the gaugino mass param- 
eter M2. We fix the top quark mass at mt = 175 
GeV. We note that this definition of the MSSM 
models is the same as in 

We adopt a phenomenological approach and al- 
low for general variations of parameters in the 
MSSM, still of course consistent with experimen- 
tal bounds and giving correct low-energy symme- 
try breaking. We keep only models that satisfy 
accelerator constraints, including the 95% C.L. 
limits on the b — + S7 branching ratio from the 
CLEO experiment §, 1.0 • lO"'' < BR(b ^ 37) 
< 3.4 • 10-"^. We calculate BR(b ^ 37) accord- 
ing to 1^. For a general supersymmetric model, 
b — > S7 receives contributions from W-bosons, 
charged Higgs bosons, charginos, gluinos and neu- 
tralinos. For our no-FCNC ansatz, the latter two 
are absent. We also impose the cosmological con- 
straint ri^/i^ < 1, where the relic neutralino den- 
sity Vl^h'^ is calculated as in pc| ]. 

We compute direct detection rates, integrated 
over deposited energy with no energy threshold, 
for pure germanium (^^Ge), sapphire (AI2O3) and 
sodium iodide (Nal) detectors. For the local 
galactic neutralino velocity distribution we as- 
sume a truncated gaussian of velocity dispersion 
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120 km/s. We fix the relative Earth-halo speed 
at 264 km/s (a yearly average). We adopt a lo- 
cal dark matter density of 0.3 GeV/cm^ when- 
ever ri-^ > f2gai, the minimum value for which 
neutralinos could make up the totality of galac- 
tic dark matter. When fi^ < Ogai, we scale the 
local density proportionally to ily^/ilsai- We use 



a heavy-squark effective lagrangian for the 
neutralino-nucleon interaction, and we use sim- 



ple gaussian nuclear form factors |12|, which are 
quite adequate for our purposes. 

We produce two initial samples of 4000 points 
each, generating model parameters randomly 
within the following bounds: ^ [100,3000] 
GeV, Ab.,At G [-3mf,3mj], tan/? e [1.2,50], 
iTiA < 1000 GeV, and, for the first sample, 
/i,M2 e [-5000,5000] GeV, while for the sec- 
ond sample we trade /Lt, M2 with the lightest neu- 
trahno mass G [—5000, 5000] GeV and its 
gaugino fraction e [0.00001,0.99999]. In these 
samples, we find no points with interesting detec- 
tion rates, say above 1 event/kg/day in Ge. But 
we notice that the distributions of points obtained 
depends crucially on our choice of sampling. So 
we do not conclude that in our class of models the 
b ^ S7 constraint is so strong to exclude inter- 
esting detection rates. We instead perform two 
special scans aiming at large neutralino-proton 
cross sections. In the first, we sample in ^-A'h as 
before but restrict niA S [0, 60] GeV. In the sec- 
ond, we sample in m^-Zg with the restricted uia 
range and further demand G [800, 1200] GeV 
and Zg e [0.01,0.99]. The results of these spe- 
cial scans for Ge, AI2O3 and Nal are shown in 
fig. |l|, together with (the top parts of) the ini- 
tial naive samplings. Remarkably, the high-rate 
zones, empty in the initial scans, are now filled 
with points. Particularly striking is the concen- 
tration around rriy^ — 1000 GeV, which obviously 
comes from the second special sampling. 

Therefore, contrary to [^, we find models that 
do not violate the experimental and cosmological 
constraints mentioned above and in which, thanks 
to a light neutral Higgs boson, the integrated 
counting rate is as large as 10 events/kg/day in 
Ge and even 100 events/kg/day in Nal. These 
models have a relatively light (100-200 GeV) 
charged Higgs boson, but they are compatible 
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Figure 1. Integrated direct detection rate R in 
^^Ge, AI2O3 and Nal detectors versus neutralino 
mass TO-v 



with the b — » S7 constraint because the charged 
Higgs contribution to the b — > S7 amplitude is 
effectively canceled, at large tan /3 and large top 
squark mixing, by the contribution from a light 
charged higgsino (cfr. ^]). In terms of the super- 
symmetry parameters this occurs when IM2I ^ 
ImI ^ ''^w and n{At + /zcot/3) < 0. We can sat- 
isfy these conditions for both positive and nega- 
tive values of fi, because we are free to choose the 
sign and magnitude of At. This freedom would 
be lost in models imposing additional theoretical 
constraints, for example in no-scale models or in 
models with a flat Kahler manifold {A = 01 
A = B — m respectively at the unification scale) . 

In our simple-minded prescription for the local 
galactic neutralino density, the value of f2gai is 
quite uncertain, both because of uncertainties in 
the density and extension of galactic halos and 
because of the poorly known relation between 
the universally-averaged and the local dark mat- 
ter densities. Fig. |2| shows the predicted rates in 
Ge versus the calculated neutralino relic density 
il^h^. The two initial samples show at the bot- 
tom and the two special samples are the band and 
cloud in the upper parts. Models with Q^h"^ > 1 
are plotted only to illustrate the trend of R ver- 
sus n^h^. Rates to the left of the vertical dashed 
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Figure 2. Integrated direct detection rate R off 
^^Ge versus neutralino relic density ^x^^- 



line {fix = ^gai) are ^-suppressed. We have cho- 
sen figai/i^ = 0.025, as in but the actual value 
might be even one order of magnitude larger or 
smaller. Were it smaller, many models in which 
the detection rate is suppressed just because fi^/i^ 
is too small would have important detection rates. 
Notice that in some of the interesting models, the 
neutralino relic density is large enough for them 
to fill up galactic halos. 

We stress that the density of points in the plots 
depends on the a priori distribution of the model 
parameters, which is entirely at our choice. One 
should not be misled in thinking that rate values 
in zones where there are more points are more 
probable than those in which there are few. This 
is evident for our special scans, for which high 
detection rates with acceptable b — > S7 branching 
ratios look "generic." No probability should be 
attached to the plotted point distributions or to 
histograms derived from them: the figures can 
only illustrate possible neutralino detection rates. 

The conservative approach we propose is to re- 
gard the whole range of outcomes of a calculation 
as a priori equally probable, irrespective of the 
parametrization. This means that really only up- 
per and lower limits can be given. Unfortunately, 
a thorough and fine scanning of parameter space 
is computationally very expensive and an alterna- 



tive analytical extremization seems prohibitively 
complicated. We therefore have to leave the fol- 
lowing question open: are there in fact additional, 
allowed points in the empty regions of our plots? 
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